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The high-mass end of the halo mass function is a sensitive probe of primordial non-Gaussianity 
(NG). In a recent study [j| we have computed the NG halo mass function in the context of the 
Excursion Set theory and shown that the primordial NG imprint is coupled to that induced by 
the non- linear collapse of dark matter halos. We also found an excellent agreement with N-body 
simulation results. Here, we perform a more accurate computation which accounts for the interval 
validity of the bispectrum expansion to next-to-leading order and extend the calculation to the case 
of a non-vanishing primordial trispectrum. 



I. INTRODUCTION 

The Excursion Set Theory initially introduced by Bond 
et al. [l| provides a self-consistent mathematical frame- 
work to infer the properties of the halo mass distribution 
from the statistics of the initial density field. The for- 
malism generalizes the original Press-Schechter idea @ 
by formulating the halo mass counting problem as one of 
stochastic calculus. The starting point is the realization 
that at any location in space the linear matter density 
fluctuation field performs a random walk as function of a 
filtering scale R. In average, this scale naturally defines 
a mass scale M — pV(R), where p is the mean mat- 
ter density and V(R) the enclosed spatial volume. By 
counting the number of trajectories which first-cross a 
collapse threshold it is then possible to compute the frac- 
tion of mass elements in halos F(M) and consequently 
derive the halo mass function dn/dM = (1/V) dF/dM, 
The requirement of first-crossing is key to solving the 
so called "cloud-in-cloud" problem affecting the original 
Press-Schechter result. In fact, the first-crossing condi- 
tion guarantees that in the small scale limit (R — > 0) and 
independently of the properties of the random walks the 
fraction of mass into collapsed objects always tends to 
unity. 

The Excursion Set is analytically solvable in the case 
of uncorrelated (Markov) random walks for which the 
evaluation of the first-crossing distribution is reduced to 
solving a standard diffusion problem. However, uncorre- 
lated random walks are generated by a special filtering 
of the linear density field which corresponds to a non- 
physical halo mass definition. In contrast any filtering 
which specifies a physically meaningful mass generates 
correlated random walks for which F(M) can be inferred 
only through a numerical computation. This has repre- 
sented a major limitation since Monte Carlo simulations 
are computationally expensive and moreover do not pro- 
vide the same level of physical insight of analytic solu- 
tions. The seminal work by Maggiore & Riotto @ has 
made a major step forward in this direction. Using the 



path-integral formulation of the Excursion Set Theory 
the authors have shown that the first-crossing distribu- 
tion of correlated random walks can be computed as a 
perturbative expansion about the Markovian solution. 
Using this methodology it has been possible to derive 
analytical formulae for the halo mass function under dif- 
ferent halo collapse model assumptions as well as Gaus- 
sian and non-Gaussian (NG) initial conditions 

In a series of papers @-IH we have used this formalism 
to evaluate the imprint of the non-spherical collapse of 
halos on the mass function. To this end we have intro- 
duced an effective stochastic Diffusive Drifting Barrier 
(DDB) model which parametrizes the main features of 
the ellipsoidal collapse of halos. Accounting for such ef- 
fects can reproduce the halo mass function from N-body 
simulations with remarkable accuracy both for Gaussian 
and non-Gaussian initial conditions. 

Here, we extend the work presented Q to derive a more 
accurate expression of the contribution to the halo mass 
function of the primordial bispectrum expanded in the 
large scale limit to next-to-leading order and compute the 
leading order contribution of the primordial trispectrum. 

The paper is organized as follows. In Section [TT] we re- 
view the path-integral formulation of the Excursion Set 
and its application to Gaussian and non-Gaussian initial 
conditions in the case of the DDB model. In Section IIIII 
we evaluate a lower limit on the interval validity of the 
bispectrum expansion at next-to-leading order. In Sec- 
tion |IV] and [V] compute the bispectrum and trispectrum 
contribution to the mass function respectively. Finally, 
we present our conclusion in Section IVI1 



II. 



EXCURSION SET MASS FUNCTION AND 
DIFFUSIVE DRIFTING BARRIER 



Here, we briefly review the main features of the path- 
integral formulation of the Excursion Set Theory. First, 
let us introduce the variance of the linear density field <5 
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filtered on a scale R: 



a 2 (R) = S(R) 



2tt 2 



dkk 2 P(k)W 2 (k,R), (1) 



where P{k) is the linear matter power spectrum and 
W 2 (k,R) is the Fourier transform of the filter function 
in real space, W(x,R). As already mentioned, by select- 
ing a volume V(R) — J d 3 xW(x,R) the filter naturally 
associates a mass to the enclosed region, M — pV(R). 
Thus we have a one-to-one relation between R (or M) 
and 5*. 

In the Excursion Set Theory the filtered density field 
S(x, R) at any random point in space performs a random 
walk as function of R, M or equivalently S which plays 
the role of a pseudo-time variable. The random walks 
start at S = with <5 = 0, since in the large scale limit 
(R — > oo) we have S — > and the matter density dis- 
tribution tends toward homogeneity, i.e. 5—^0. We are 
interested in counting trajectories that at a given value 
of S cross for the first time a collapse density threshold 
B such that 5 = B. This threshold encodes all informa- 
tions on the gravitational collapse of halos. In order to 
model features of the non-spherical collapse, the absorb- 
ing barrier B is promoted to a stochastic variable (see 
e.g. 0, E3, E3) also performing a random walk as func- 
tion of S (i.e. R or M). In such a case it is convenient 
to introduce Y = B — 5, which performs a random walk 
starting at Y(0) = Yq with barrier crossing at Y c = 0. 

Our goal is to compute the probability distribution 
n(Yo, Y, S) of trajectories starting at Yq which reach the 
value Y at S without ever touching the barrier Y c = 0. 
This can be computed as a path-integral over the en- 
semble of the random trajectories (see Q for a detailed 
derivation). Let us discretize the pseudo-time variable S 
in equally spaced steps, AS = e, such that Sk = ke with 
k = l,..,n. The probability distribution of trajectories 
starting at Yo and ending in Y n at S n and that have never 
crossed the barrier before is given by 



U £ (Y a ,Y n ,S n )= dYi.. dY n ^W(Y ,..,Y n ,S n ), 

(2) 

where 



W(Y ,..,Y n ,S n ) = / V\e^=^ Y *{e £1 



where the brackets (...) refer to an ensemble average of 
the random walks and the averaged quantity is the ex- 
ponential of 

z = E E - E A - -K to -Y ip ) c (4) 

p=l ii = l i p =l 

which is the partition function of the system written 
in terms of the p-point connected correlation functions 
(Yi 1 ..Yi p ) c of the random walks. Thus, the the proper- 
ties of the stochastic system are entirely determined by 



the connected correlators. Once these are specified then 
Eq. ^ can be integrated in the continuous limit to finally 
obtain the first-crossing distribution 



/OO 
dYU(Y a ,Y,S) 



and the halo mass function is given by 

dn _ p dlogCT -1 

dM ~ ■'^'JP dlogM ' 



(5) 



(6) 



where f(a) = 2SJ-(S) is the so called "multiplicity" func- 
tion. 



A. Gaussian Initial Conditions 

Let us consider a Gaussian density field smoothed with 
a top-hat filter in real space. On average the density 
field is homogeneous, thus implying that {6(S)) C = 0. 
Therefore, due to the Gaussian nature of the field the 
only non-vanishing connected correlator is the 2-point 
function (S(S)5(S')) C , while all higher-order connected 
correlators identically vanish. Maggiore & Riotto [3] 
have shown that for standard cosmological scenarios 
with Cold Dark Matter power spectra, the 2-point func- 
tion smoothed with a sharp-x is well approximated by 
(S(S)S(S')) C = mm(S,S') + A(S,S'), where the first 
term corresponds to Markov random walks generated by 
a sharp-k filter and the second term is well approximated 
by A(5, S') = kS/S'(S' — S) with a nearly constant am- 
plitude k < 1. Thus, the pseudo-time correlations in- 
duced by the filter function can be treated as small cor- 
rection about the Markovian case and the mass function 
obtained using a perturbative expansion of the partition 
function in the path-integral in powers of k. 

Concerning the barrier random walks, in we have 
introduced a stochastic model with linear drift and Gaus- 
sian diffusion characterized by (B(S)) = 5 C + j3S and 
(B(S)B(S')) C = D B mm(S,S'), where S c is the linearly 
extrapolated critical spherical collapse density, j3 is the 
average linear rate of deviation from the spherical col- 
lapse prediction and Db is the amplitude of the scatter 
about the average . In such a case the non- vanishing 



1 In the Excursion Set the barrier diffusion coefficient parametrizes 
the stochasticity inherent to the ellipsoidal collapse of halos. 
However, it is important to keep in mind that in the Excursion 
Set halos can form out of any random position. On the other 
hand, numerical simulations show that halos form preferentially 
out of peaks of the linear density field as suggested by the hi- 
erarchical model of structure formation. Thus when comparing 
with N-body results the value of Dg can be biased by the un- 
derlying assumption of the Excursion Set approach (see ll'J for 
an extension of the formalism to random walks around density 
peaks). 
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connected correlators of the Y variable are 



-(6(Si)S(Sj)S(S k )) c with 



(Y(S)) C = S C + /3S (7) 
(Y(S)Y(S')) C = (l + D B )min(S,S')+A{S,S').{8) 

Substituting these expressions in Eq. Q and perform- 
ing a double expansion in n and /? we have derived the 
Gaussian multiplicity function 



f G (a) = f (a) + f K=1 (a) 



(9) 



where /o(c) is the Markovian contribution and / K= i(cr) 
is the filter correction to first order in n and up to second 
order in j3 which read as 



S c 2a _ 
Jo(cr) = — \ — e 
a V 7r 



(10) 



and 



/ K= i(<7) = /"(a) + fZ^{a) + f™-£(a) (11) 



with 

rm — m ( \ 



-K — 

a 



e 2« 



Wo, 4 

2 V 2d 2 



/l"(^)=-«^/? 



ftErfc £ c 



7x^(2) (o-J = -a/3 



ft J2rm—rnr_\ . r rm — m / \ 



(13) 



(14) 



where k = ok and a = 1/(1 + Dg). Equation (TI4")) 
includes a term C(/3 2 ) which was missing in the orig- 
inal derivation presented in 0, As we explain in 
Appendix [A] this is due to having neglected a fac- 
tor exp(— (3 2 S/2) in the computation of the probabil- 
ity n e (Y c , Y C ,S) 2 and which enters the calculation of the 
memory-of-memory term to first order in k. 



(*(5 i )*(S i )*(S fc )>c = 

x W(kj,Rj[Sj])W(k 
x (C(k i )C(k i )C(kfc)>c 



/c^ d k-j d kfc r ~ 



W(ki,Ri[Si]) 



»"^»" B I (2tt) 3 (2tt) 3 (2tt) 3 
x IF(fc, , i?, [5,-] ) VK(fc fc , i? fc [S k ])M (h)M fa )M (kj)x 



(15) 



where W is the Fourier transform of the sharp-x filter, 
M{k) = 2/(5H$n m )T(k)k 2 , H is the Hubble constant, 
VL m the matter density, T(fc) the transfer function and 
C(k) is the curvature perturbation with 

(C(k l )C(k,)C(k fc )>c = (2n) 3 5 D (k l + k J +k k )B(h,k,,k k ), 

(16) 

where B(ki, kj, k k ) is the so called "reduced" bispectrum. 

By expanding Eq. ([3]) in powers of the amplitude of 
the reduced bispectrum (usually parametrized by the co- 
efficient /jvl) ; we obtain to first-order in Jnl the non- 
Gaussian part of the first-crossing distribution 



d 

Fng{S) = —-ggF NG (S) 
where Fng(S) is the continuous limit of 



(17) 



n 

( 12 ) F NG (S) = - ]T (6(Si)5(Sj)8(S k ))t 



i,j,k=Q 

fOO poo 



(18) 



x dY dY 1 ...dY n - 1 d i d j d k W Q {Y ,...,Y,S), 

JY a JY C 

where Wq(...) is the Gaussian Markovian probability den- 
sity distribution. Eq. (1181) can be evaluated provided we 
have an analytical expression for the primordial bispec- 
trum. In @ we have used the standard approach of con- 
sidering a triple Taylor series of the primordial bispec- 
trum in the large scale limit. In the next Section we will 
study in detail the range of validity of such an expansion 
and infer the relevant contribution to the non-Gaussian 
halo mass function. 



III. INTERVAL VALIDITY OF PRIMORDIAL 
BISPECTRUM EXPANSION 



B. Non-Gaussian Initial Conditions 

In the case of non-Gaussian initial conditions the 
higher-order connect correlators of the linear density 
field are non- vanishing. Let us consider the case of pri- 
mordial non-Gaussianity sourced by a bispectrum term, 
hence in addition to Eq. (JT]) and (JU), the partition func- 
tion contains the contribution of a non- vanishing 3-point 
connected correlation function (Y(Si)Y(Sj)Y(S k )) c = 



Let us expand the bispectrum Eq. (|15[) in a triple Tay- 
lor series around Si ~ Sj ~ S k ~ S (for convenience we 
set S = S n ): 



plq\r\ 



p^q^r—O 

x (S-SiY(S-S k ) r C^' 9 ' r) (S) 



where 

G^ r \s) ee ^^^-(SiS^iSMSk)), 



We thank Ruben van Drongelen for pointing this to us. 



(19) 



i,j,k=n 

(20) 
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We expect the signature of primordial non-Gaussianity 
to be stronger at large scales (S — > 0) where the evolution 
of the density field remains linear. Hence, the leading 
order contribution to the halo mass function is given by 
the lowest order term of the bispectrum expansion. This 
corresponds to having p+q+r = which gives the leading 
order term (S 3 (S)) that can be computed numerically 
using Eq. (|15[) for a given type of primordial NG. It is 
convenient to introduce the normalized skewness S 3 (S) = 
(S 3 (S)}/S 2 . In Q we have provided fitting formula for 
^3(5) and its derivatives accurate to a few percent for 
local and equilateral non-Gaussianity. 

The next-to-leading order contribution is given by 
three terms corresponding to the case p + q + r = 1. 
Hence, up to next-to-leading order the primordial bis- 
pectrum reads as 

(S(Si)5(Sj)5(S k )) c = (S 3 (S)) - (S - S t ) G { i'°-°\S)+ 



-(S-S 3 )G^ 0) 



(S)-(S-S k )G^°' 1) (S), 



(21) 



since Si ~ Sj ~ Sfc we can collect the terms in (S — Si), 



moreover by computing G 3 ' , G 



(1,0,0) .-,(0,1,0) 



and G 



(0,0,1) 



as 



derivatives of Eq. ([15)) one can notice that 

dR d 



Gl' ' (S) + G /' (S) + G /'\S) 



dS dR 



and introducing 



(S 3 (R)) (22) 



(23) 



we can rewrite Eq. (|21[) as 

(*(5i)*(5 i )*(S fc )>c « S 2 S 3 (S) - (S - Si)SU 3 (S). (24) 



We can now derive a lower limit, S n 



on the value 



of Si,Sj,Sk for which such an expansion remains valid. 
This is obtained by imposing the next-to-leading order 
term to be smaller than the leading one. We find 



Smin — S 



1 - 



U 3 {S) 



(25) 



using the fitting formulae for S 3 (S) and U 3 (S) derived in 
9] we find that to good approximation S m i n = S/a with 
a^ c = 0.373 and aC ui = 0.382 for local and equilateral 
NG respectively. 



IV. NON-GAUSSIAN HALO MASS FUNCTION 
AND BISPECTRUM EXPANSION ACCURACY 

Having inferred a lower limit on the interval validity 
of the bispectrum expansion up to next-to-leading order 
we can infer a more accurate estimate of its contribu- 
tion to the multiplicity function. Given the bispectrum 
expansion Eq. ([24]) we can split Eq. (fT8|) as 

Fng(S) = Fk G (S) + F^(S). (26) 
As shown in [9| the leading order term is given by 



1 f°° d 3 

f^ g (s) = -s 2 s 3 (s) J y dY_n (Y ,y,s) 



(27) 



where Uo(Yo,Y, S) is the probability distribution of the 
Gaussian random walks in the case of the Diffusing Drift- 
ing Barrier model given by Eq. (8) in Q. The integral 
can be computed analytically to finally obtain the leading 
order contribution to the multiplicity function 



a 1 2a 

D V 7T 



a(S c + fia 



2\2 



2a 2 



cr 4 cr <y 



dS 3 (a) 
dloga 



—S 2 - 1 + 3a/3£ c + 4aP 2 a 2 



-a 3 /3 3 cr 4 e - 2a/3<5 <=Erfc 
3 



^(Sc-f3a 2 ) 



453(a) + 



dS 3 (a) 
dloger 
(28) 



On the other hand let us detail more the derivation of Q . In such a case we have from Eq. ([T8|) that 
the next-to-leading order term which differs from that of 



F^(S) = -±SU 3 (S)J2(S-S i )x 

/CO poo poo 

dYY, dYt... dYn-ididjdkWo, 
j,k 

(29) 



5 



where we have decomposed the sum in Eq. (1181) and kept 
only the terms up to n — 1. In fact, as shown in [9( 
the integral in dY n of the terms with i,j, k = n vanishes 
since the integrands are total derivatives. Furthermore it 
is easy to show that Ylj k d 2 /dY 2 , thus 



n — l 

Fn8(s) = -^su 3 (s) J2( s ~ s * 



dY 2 



dY x . 



dY^diWo 



(30) 



and 



g d 2 lNL 

dS dY? 



dS * c -^r(y°+^i) 2 x 



s 



3/2 



x <^ -2aj3 + a 2 /3 2 S c + a 



Si 



3 + 2af36 c + a 



Si 



e 2 



2V5 - SS 



l-I + I- * 



a a S — Si 



1 q2 5 + 5, 

—an — , 

a H y/S^Sl 



2 V a 



-0 x l-(S-Si) 



(35) 



where the sum is bounded from below due to the fact 
that S m i n < Si < S. The multiple integral in the above 
expression can be computed by part and using the fact 
that Wq obeys the Chapman-Kolmogorov equation 



Wo(Y ,..,Y i ,..,Y,S)=Wo(Y ,..,Y c ,S i -S itain )x 
x W (Y C ,..,Y,S-Si + S imin ) 

we have 

1 n— 1 

F$8 = -xSUs(S) J2 (S-Si)x 



(31) 



d 2 

< ly 



n ( y , y c , Si - - )n(V C) y, s - s, + - 

a J \ a 
(32) 



-2x10" 



-4x10- 




local NG 



0.3 



0.4 



where U(Y ,Y C , Si ~ S/a) and IL(Y C , Y.S— Si + S/a) arc 
given by Eq. (3.20) and (3.21) in Q respectively. In 
the continuous limit and taking into account that the 
bispectrum expansion up to next-to-leading order is valid 
in the range S m i n ~ S/a < Si < S we have X^2=i ~~ * 

Is /a dS*- 

Finally, we obtain the first-crossing distribution to 
next-to-leading order which reads as 



J~ng(S) 
with 



° 1 NL 





\d 2 I% L 


u 3 - 




6 ' 


I dY 2 







qu d d 2 I% L 
SU3 dS^Yf 
(33) 



dY 2 



dS„ 



ft J s. s 



3/2 









-Si 







-3 + 2a/38 c + a-£- 

Si 



-2a/3 + a 2 p 2 5 c + a-£ 

Si 

, , , r ffl a (S-S0 + _ 5i )Erf c 



(34) 



5xl0" 6 



-5xl0- 8 




equilateral NG 



0.3 

lOg 10 (<T- 1 ) 



0.4 



FIG. 1: Leading order (red dash line), next-to- leading order 
(blue dot line) contribution to the non-Gaussian multiplicity 
function (black solid line) for local (top panel) and equilateral 
(bottom panel) non-Gaussianities respectively. 

The integrals in Eq. and (|3"5j) can be computed 
numerically to finally evaluate the next-to-leading order 
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contribution to the multiplicity function, 



V. TRISPECTRUM CONTRIBUTION TO THE 
NON-GAUSSIAN HALO MASS FUNCTION 



JnGV*) = 2(1 S- NG {<J ) 



(36) 



In Fig. (Q]wc plot /jvg( <7 )' /jvg( (T ) an0 - their sum in 
units of Jnl for local (top panel) and equilateral (bottom 
panel) non-Gaussianity. We have set the DDB model pa- 
rameters to the values best fitting the Gaussian halo mass 
function inferred from Gaussian N-body simulations in 
[l3j . The mass interval shown here corresponds to that 
probed by these numerical simulations. We may notice 
that in the local and equilateral cases the leading or- 
der term is larger than the next-to-leading order one at 
high-masses (log 10 cr _1 > 0.1), while at lower masses the 
next-to-leading order is larger. This is expected since 
as already mentioned the signature of primordial non- 
Gaussianity at large masses results of the lowest order in 
the bispectrum expansion. 

We can now evaluate the overall contribution to the 
halo multiplicity function, f(a) — /g(c) + /vg(o')- In 
Fig.[5]we plot the relative difference of the NG halo mass 
function with and without next-to-leading order term for 
local (top panel) and (bottom panel) equilateral non- 
Gaussianity respectively in the case of Jnl = 150. As 
we can see the differences is no larger than 2% in the low 
mass range, hence the next-to-leading term remains neg- 
ligible even for large non-Gaussianities in the mass range 
corresponding to halos with M > 1O 13 M0 and can be 
neglected for practical purposes. 




A number of scenarios of primordial inflation predict 
deviation from Gaussianity of the form (see e.g. flil Il5j|) 

C = Cg + lfNL(C 2 G - (Cg» + 7^9nl( 3 g + 0(C 4 G ), (37) 

where (g is the primordial curvature perturbation and 
9nl is the amplitude of the cubic term which give rises 
to a non- vanishing 4-point connected correlation function 
of the linear density field given by 

W(h, Ri)W[k h Rj)W(k k , R k )W{k u R t )x 

x^(fc i )M(%)7W(fc fe )A4(fcO(C(k < )C(k J )C(kfe)C(kO)c 

(38) 

where (C(k l )C(k,)C(k fc )C(k i ))c = (27r) 3 fe(k 4 + k, + k k + 
ki)T(ki, kj,kk,ki) and T(ki, kj 7 kk, h) is the trispectrum. 
The amplitude of the trispectrum is affected not only by 
the cubic term in Eq. (|37[) and parametrized in terms 
of gNL, but also from the skewness as parametrized by 
/jvl- In particular, For models in which curvature per- 
turbations are sourced by a single scalar field, as in the 
case of the curvaton model (see e.g. 16]) the skewness 
contribution to the kurtois is given by rff L — 36/25/^ L . 

We can compute the trispectrum contribution to the 
multiplicity function by including the 4-point connect 
correlator in the partition function and expand the path- 
integral for small values of the trispectrum amplitude. As 
in the case of the bispectrum, to leading order in a large 
scale expansion the trispectrum can be approximated as 

${SiWiW*Wi))* - ^\S)) C (39) 
To first-order in the trispectrum amplitude we have 



(40) 



where F^g L (S') is the continuous limit of 



i,j,k,l 



dY J dY 1 ...dY n -id i d j d k diW . 



(41) 



FIG. 2: Relative difference bewteen of the non-Gaussian halo 
mass function with and without next-leading order contri- 
bution in the case of local (panel a) and equilateral non- 
Gaussianity (panel b). 



Using the fact that ■ k l — > d 4 /dY* we obtain 



Il (Y ,Y,S)dY, 



c JY, 



(42) 
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FIG. 3: Leading order contribution of the bispectrum (red 
dash line), the trispectrum (blue dot line) and their sum for 
local type of primordial non-Gaussianity. 



the integral can be computed analytically, 

_l 6(a/3) 3 + 8 ^ + q^ Yq _ U «L(P 5c - 8^(36* + 



2-5* 
S 3 



i-2a/3Y 



Erfc 



(43) 



As in the case of the bispectrum, it is convenient to 
introduce the reduced cumulant of order 4 



s 3 ' 



(44) 



substituting in Eq. (l4"2j) and evaluating the first-crossing 
distribution Eq. f|40[) we finally obtain the trispectrum 



contribution to the multiplicity function: 



nj4 -2a/35 c 



x Erfc 



x (acr) 2 



3> -/*>*) 





3 dloger 
dloga 



r(8 c +f3a 2 ) 2 . 



6 3 8 



(45) 



It can be noticed that by setting the barrier model 
parameters to the spherical collapse values a = 1 and 
ft = we recover the formula derived in [17] . It is also 
worth noticing that in the spherical collapse limit and 
neglecting the filter correction to first order in k, the NG 
multiplicity function given by the sum of the Markovian 
term, the bispectrum and trispectrum leading order con- 
tributions has the same functional form as that derived 
in [l8| using the Edgeworth expansion to describe the 
non-Gaussian probability distribution of the initial den- 
sity perturbations. 

In Fig. [3] we plot the multiplicity function for the bis- 
pectrum and trispectrum respectively for different values 
°f Inl and c/jvl within Cosmic Microwave Background 
limits (20[, having used the the fitting formulae of 54(17) 
derived in and set the DDB model parameters to 
their Gaussian value, which is a good approximation for 
/nl < 150. We can see that even for g NL = 10 6 and 
tnl = 10 4 the trispectrum exceeds the bispectrum sig- 
nal only in a very limited low mass range. In contrast 
for gNL — and tnl — 10 4 the effect on the non- 
Gaussian multiplicity function remains very small com- 
pared to that of the bispectrum. 



VI. CONCLUSION 

The halo mass function carries an imprint of the statis- 
tics of the primordial density field as well as the proper- 
ties of the halo collapse process. The path-integral for- 
mulation of the Excursion Set theory provides a pow- 
erful and self-consistent mathematical framework to ac- 
count for these effects on the halo mass function. Here, 
we have extended a previous analysis [i| and performed 
a more accurate derivation of the contribution of the 
primordial bispectrum expanded in the large scale limit 
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to next-to-leading order for the Diffusive Drifting Bar- 
rier model introduced in 0, H|. We have shown that 
the next-to-leading order term of the primordial bispec- 
trum decomposition contributes to no more than ~ 2% 
of the non-Gaussian mass function. Thus, for all practi- 
cal purposes it can be neglected. We have also derived 
an analytic formula for the trispectrum contribution. As 
in the case of the bispectrum, the multiplicity function 
depends on terms which couple the parameters encod- 
ing the ellipsoidal collapse of halos with the primordial 
four-point correlation function. Also in this case we find 
that in the spherical collapse limit the trispectrum con- 
tribution reduces to the functional form derived in the 
Press-Schechter formalism using the Edgeworth expan- 
sion. However, in order to reproduce N-body simulation 
results the latter requires two ad-hoc prescriptions. First, 
the non-Gaussian prediction is rescaled by a Gaussian 
simulation calibrated multiplicity function, such as to ac- 
count for the imprints of the ellipsoidal collapse. Thus, 
implicitly assuming that the effect of the non-spherical 
collapse of halos on the mass function is independent of 
the amplitude of primordial non-Gaussianity. However, 
this is not the case since our derivation explicitly show 
a coupling between the halo collapse parameters and the 
non-Gaussian contribution. Then, an empirical "fudge 
factor" is introduced to correctly reproduce the shape 
of non-Gaussian mass function from N-body simulations. 
Hence, it is not surprising that such analyses find that 
for large primordial non-Gaussianities the fudge factor 
depends on /jv l ■ The Excursion Set derivation provided 
here shows that these two empirical prescriptions are not 
necessary and provides a better physical insight on how 
the stastics of the initial density field and the non-linear 
collapse of halos influence the mass function over the 
range of mass which will be probed by the upcoming 
generation of cluster surveys. 



ing at the barrier, H £ (Y C , Y c , S n ) which enters the calcu- 
lation of the memory-of-memory term due to the non- 
Mar kovian filter corrections (see Eg. (11) in !7|). Fol- 
lowing the derivation presented in 8] let us consider the 
Chapman-Kolmogorov equation: 



x n e (Y"o> Yn-1, S n -i), 



(Al) 



with 



(A2) 



expanding in Taylor series the left-hand-side of Eq. (IA1|) 
in powers of e and the right-hand-side in powers of Ay 
we obtain 



: ds: 



2 dsi 



+ ... = 



dx\ — 

V a 



Bp 



(A3) 



It is worth noticing that the Markovian probability dis- 
tribution starting at the barrier value Yq = Yc vanishes 
in y/e while for Y = Yq = Y c the first finite corrections 
to the probability distribution function is in s. At order 
one in e for Y = Y c , this equation gives 



dS n 



ad 2 IL F 



2 dY 2 



pan, 

2~dY ' 



(A4) 
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Appendix A 

Here, we derive the finite correction to the probabil- 
ity distribution of the random walks starting and end- 



which has the same form of than the Fokkcr-Planck equa- 
tion associated with the Markovian solution for the DDB 
model in the continuous limit whose solution for a barrier 
at a generic point Y c is given by Eq. (3.11) in In such 
case for Y = Yq = Y C the drifting term stay in factor of 
Gaussian minus anti-Gaussian. Therefore we can assume 
that 



U e (Y c ,Y c ,S n ) = Ce- 



-a/3 2 S n /2 



S. 



3/2 



(A5) 



The value of the constant C can be evaluating using 
the path-integral Eq. ([2]) for n = 2 with Yq — Y^ = Y c and 
the explicit form of Wq(..) given by Eq. (A3) in Thus, 
equating Eq. (|A5[) on the left-hand-side to the result of 
the integral over Wq(...) on the right-hand-side we obtain: 

C= x/^. 
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